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ABSTRACT 


A propagation model which combines, a Gaussian profile for particle 
release from the sun, with interplanetary particle densities predicted 
by focused diffusion, is proposed to explain the propagation history of 
flare-associated energetic particles. This model, which depends on 
only two parameters, successfully describes the time -intensity profiles 
of 30 proton and electron events originating from the western hemisphere 
of the sun. On comparing the predicted and observed density profiles, 
it is found that the scattering mean free path X is on the order of 0.1 
to 0.3 AU for 4-80 MeV protons and is a factor of 2 to 3 smaller for 
0. 5-1.1 MeV and 3-12 MeV electrons. Since the rigidities here differ by 
more than two orders of magnitude, this implies that X is only weakly 
dependent on rigidity. The rigidity dependence of X, as well as its 
magnitude, strongly disagrees with that predicted by the quasilinear theory 
of pitch-angle scattering. Generally, particles are released from the sun 
over a finite interval. Approximating this profile of injection by a 
truncated Gaussian, the r.m.s. width a is less than 1 hour for most events. 
The injection width a decreases with velocity, but it does not depend 
on rigidity. In almost all events, particle release begins at the time 
of flare acceleration. 

Subject headings: cosmic rays: general - hydromagnetics - 

interplanetary medium - sun : flares 
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I. INTRODUCTION 

Although the behavior of charged particles in interplanetary space 
has been studied for more than two decades, the exact nature of their prop- 
agation is not yet clearly understood. Classically, the scattering of 
charged particles in random magnetic fields was treated in terms of dif- 
fusion (Parker 1956, 1963; Meyer, Parker, and Simpson 1956; Bryant et al . 

1962; Jokipii 1966, 1967; Burlaga 1967; RoeloJ 1966; 196Sl Englade 1971; Lupton 
and Stone 1973,- Webb et al. 1973; Webb and Quenby 1973). Recent work on solar- 
particle events has reached a consensus that the scattering mean free path of 
low energy protons is on the order of 0.1-0. 2 AIT (Lanzerotti et ad. 1973; Roelof 
and Krimigis 1973 r Reinhard and Wibberenz 197*+; Ma Sung, Van Ho-llebeke , and 
McDonald 1975; Countee and Lanzerotti 1976; Zvickl 1976; Hamilton 1977- ; Ma Sung 
1977 ) s which is an order- of magnitude larger than the mean free path estimated 
from the interplanetary magnetic power spectrum by Jokipii and Coleman (1968) . 

In the regime of weak scattering, implied by this large mean free path, 'the 
treatment of particle propagation in terms of simple diffusion is not appro- 
priate. Recently, (Earl 1976a, b) and: Earl and Bieber (1977) develope<Ta 
theory which treats particle transport in this regime, and which gives a 
quantitative prediction of both the intensity and the anisotropy of the par- 
ticle fluxes in space hased on a Fokker-Planck description of pitch-angle 
scattering (Jokipii 1966, 1967; Roelof 1.968).. This theory also- included the 
focusing effect which results from conservation of the first adiabatic invariant 
for particles propagating along a diverging magnetic field (Roe-lof 1966). In 
this paper, the basic propagation process of flare-associated particles in 
interplanetary space is studied in terms of a detailed model which incorporates 
a coronal injection profile and the theory of focused diffusion in inter- 
planetary space. This model assumes purely longitudinal propagation uncom- 
plicated by azimuthal transport in the corona or interplanetary space. 
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Unlike previous works on diffusion whi-ch involked an escaping boundary beyond 
1. AU to account for the observed exponential decay of particle events, this model 
treats interplanetary propagation in terms of only one parameter, (V/AL) , which 
gives the relative strength of scattering to focusing, and which predicts an exponen- 
tial decay in agreement with observations. Because the model takes into account 
non-diffusive effects which are particularly significant during the initial phase 
of an event, we have been able to deduce, by comparing the observed onset and decay 
phases with the model, not only the properties of interplanetary propagation, but 
also the characteristics of coronal escape. While the existence of a finite coronal 
injection has been pointed out previously (Feit 1973; Nielsen, Pomerantz, and 
West 1975; Palmer, Palmeira, and Allum 1975; Reinhard 1975; Ng and Gleeson 1976; 
Wibberenz, Lanzerotti, and Venkatesan 1976), it is firmly established here for 
30 events, and the injection time is demonstrated to depend on particle velocity 
but not on rigidity. 

Particle data were collected by the Goddard Space Flight Center cosmi-c ray 
experiments on IMPs IV and V. Some thirty events associated with flares on the 
western hemisphere of the sun were tested against the model. 

Section II describes the particle data and flare association. Section III 
gives the propagation model. In section IV, the profiles predicted by the model 
are compared with electron and proton events, and in section V, the results are. 
discussed. This study gives important new information on a) the rigidity and 
velocity dependence of low energy particle propagation in the interplanetary 
magnetic field, b) the velocity dependence of particle release from the sun, c) 
the time of release of particles to flare acceleration. 
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II. DATA 

Electron and proton data were obtained by the Goddard Space Flight 
Center cosmic ray experiments aboard the highly eccentric earth orbiting 
satellites IMPs IV and V. They cover five years between May 1967 and 
December 1972 with short interruptions from May 1969 to June 1969 and from 
November 1971 to February 1972. 

The GSFC instrument consisted of two dE/dx vs. E. telescopes. A low 
energy detector (LED) detected protons and alpha particles between 4 and 20 
MeV/nucleon, and a medium energy detector (MED), detected protons and alpha 
particles between 20 and SO MeV and electrons between 3 and 12 MeV. The thin 
dE/dx element, MED A, was used to detect electrons in the 0.5- 1.1 MeV energy 
range. This element had a geometric factor of 60 cm 2 -ster. The axis of the 
MED was paxall to the spacecraft spin axis which was normal to the ecliptic 
plane; the axis of the LED was normal to the spin axis with a field of view of 
<20°. The satellite spin rate was approximately 23 rpm. These detector systems 
were described in detail by Bryant et al. (1962), Kinsey (1970), and Van 
Hollebeke, Wang, and McDonald (1974). Detectors on IMP IV were essentially 
identical to those on IMP V. 

, The theory described in Section III predicts not only the isotropic par- 
ticle density, but also the detailed form of the anisotropy. Obviously, it 
would be very desirable to compare these predictions with observed anisotropies. 
Unfortunately, the present data define only the isotropic component, for they 
were derived from telescope counting rates averaged over spin azimuth. 

During the five-year period, 185 solar cosmic ray events were detected. 
Their time histories were catalogued by Van Hollebeke et al. (1974). To obtain 
a reasonable resolution above the galactic cosmic ray background and to provide 
a better separation from co-rotating events, the proton flux above 20 MeV was 
required to exceed 10 _4 particles/cm 2 -sec-sr -MeV for several hours. Among 
the 185 events detected, 125 were associated with an optical flare, a soft 
X-ray (2-20 X) burst, spectral type II, III, IV radio 
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bursts or intense microwave bursts with peak flux > 500kl 0~ 22 watts nf 2 Hz _1 
at frequencies above 3000 MHz (Van Hollebeke, Ma Sung, and McDonald 1975) . 

In this paper, we consider 70 proton and electron events, which orig- 
iginated between MV 20 and MV100 , and which comprised approximately 56% 
of all identified events. The longitudes span a region of fast propagation 
and preferred connection (Reinhard and Wibberenz 1974; Van Hollebeke et 
al . 1975, Ma Sung et al. 1975), within which coronal effects on particle 
propagation are minimal. To properly compare the particle time histories 
with model prediction, the events selected for analysis were required to 
have a sufficiently clean background (i.e. no significant flux of solar 
cosmic rays from a previous event) and to have a maximum proton intensity 
above 4'Xl0^ 4 /cni 2 -sec-sr-MeV at 40 MeV. In addition, events which showed 
distinct local perturbations or which occurred at the time of perigee 
passage were excluded. These requirements reduced the sample to a total 
of 31 events, of which 30 were found to agree with the model described in 
Section III. All but 3 events, those on January 28, 1970, January 29, 1970, 
and April 22, 1971, were well associated not only with H a and X-ray 
emissions, but also ivith type II, type IV or intense microwave bursts. 
Therefore , we believe that the flare association is sound in all except 
perhaps these 3 cases. 
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III. THE INTERPLANT ARY PROPAGATION MODEL 


The model is based on the theory of focused diffusion (Earl 1976a), 
which describes the particle distribution function f(y,z,t) in terms of- 
solutions of the one dimensional Fokker-Planck equation, 


3f * „ 3f 
3t pV 3Z - 


3_ . 9f 

8li ^ 3y 


V 2. 3f 

2L O* ) ^7 


CD 


where y is the cosine of the pitch angle, z is the distance along the 
guiding field, V is the particle velocity, and L is the scale length for 
spatial variation of the guiding field, L = -B/(3B/3z). The Fokker-Planck 
coefficient for pitch-angle scattering was assumed to have the form, 

Kh) = Alvf^ 1 (1-y 2 ). (2) 

According to quasilinear theory (Jokipii 1966, 1971), the parameter q is the 
spectral index for interplanetary field fluctuations, and the parameter 
A depends on the velocity and the rigidity of the particles, and on the 
amplitude of the random fields. Because the validity of quasilinear theory 
has been questioned in recent years, it is important to realize that this 
form was adopted, in the absence of a definitive- analysis of pitch-angle 
scattering, as a working formula whose parameters conveniently characterize 
the intensity A and the anisotropy- q of the scattering. According to 
this point of view, the magnitude and rigidity dependence of these parameters 
will be regarded here as empirical values to be determined experimentally 
rather than as calculated values related to the interplanetary fluctuation 
spectrum by quasilinear theory, . In .this way, we attempt to specify the 
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dependences that must be predicted by sound theories of pitch-angle 
scattering. 

The assumption of one dimensional propagation along the mean field 

is supported by the observation of strongly field aligned anisotropies 

during the initial phase of low energy proton and electron events (McCracken, 

Rao, and Bukata 1967; McCracken, Rao, and Ness 1968; Fan et al. 1968; Allum 

et al. 1971; Rao et al. 1971; Roelof and Krimigis 1973). Other evidence for 

field aligned propagation includes the distribution of events in helio longitude 
which shows a pronounced peak between W30 and W60 (Van Hollebeke et al. 1975), 

and the distribution of the event risetimes which shows a broad minimum- 

between W20° and W90° (Barouch, Gros and Masse 1971; Datlowe 1971; Simnett 

1971; McKibben 1972; Lanzerotti 1973; Reinhard and Wibberenz 1974; Ma Sung 

et al. 1975; Van Hollebeke et al. 1975; Ma Sung 1977). 

Except for trivial special cases, exact solutions of the transport 
equation cannot be obtained in closed form. Consequently, to obtain , 
useful approximate solutions, Earl (1976a) expanded the distribution function 
in terms of the eigenfunctions of an operator that describes the combined 
effect of scattering and focusing. This representation leads to a matrix 
formulation that is completely equivalent to the -original transport equation. 
However, this formulation in terms of infinite matrices can be drastically 
truncated with no appreciable loss in -accuracy, because only the first 
few terms of the eigenfunction series are 'significant. Thus, the present 
model is based upon analytical solutions which retain only the two lowest 
order eigenfunctions. However, the phenomenon of dispersion, which involves 
two eigenfunctions of the next higher order, was also taken into account 
approximately by the convolution method described below. 
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It seems test not to repeat here formulae- for the analytic solutions, 
for they are complicated expressions that involve many parameters. Instead, 
we .emphasize the crucial point that the exact form of these solutions 
is completely specified by only two parameters, q and (V/AL) . The latter 
parameter, which is the ratio of the scattering length (V/A) to the 
focusing length L, characterizes the importance of focusing relative to 
scattering. Throughout this work, a fixed value, q = 1.5, was adopted. 

This value approximates the observed spectral index of interplanetary 
magnetic fluctuations. Similarly, the focusing length L was fixed at its 
value at 1,2 AU, L = r/2 =0.6 AU. Thus, in effect, the only free parameter 
was the scattering length (V/A) . 

The parameter (V/A) and the index q are related to the diffusion 
coefficient K by (Earl 1974b) 

Y _ (V 2 /A) 

(2-q) (4-q) * . (3) 


Because K = (1/3)VA, where A is the mean free path, it follows that. 


A = 


3 

(2-q>(4-q) 



(4) 


Hence, A increases as (V/A) increases while q is constant, or as q increases 
while (V/A) is constant. During typical western hemisphere events, in which 
(V/AL) “0.1-0. Z as obtained in this analysis, the diffusion coefficient and 
the mean free path for 40 MeV protons (V = 1.98 ATJ/hr) range from 
(0.6 to 1.2)xl0 22 cm 2 /sec and from 0.14 to 0.28 AU, respectively. 

The theoretical density profile we compared with experimental data is 
the function F 0 + {z,t}, whose specific form, given by Earl (1976a, eq. 55), 
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will not be repeated here for the reason given- above. It was evaluated 
by substituting in Earl's equation (55] the parameters defined by his 
equations (37] to (41) which were calculated as a function of (V/A] from 
the numbers that appear in his Table I for q = 1.5. The profile. F o + {z,t} 
consists of a 6 - function coherent pulse, which moves with speed V = 0.6V 
outward alon^ the guiding field, and a wake, which trails behind the 
pulse. A second component of the general solution, the function F ~{z,t}, 
which could be significant if particles were injected inward, was not 
included because flare injection clearly corresponds to the outward injec- 
tion described by F + {z,t}. 

In a magnetic field of the form F “l/r^, the focusing length L is not 
constant, for it varies as L = r/2. Thus, as L increases with distance 
from the sun, the focusing parameter (V/A) decreases from the supercoherent 
regime (V/A)» 1, within which solar particles propagate in a well collimated 
beam, to the regime of focused diffusion V/AL <£.1, within which the present 
formulation applies. Because this transition- fronr supercoherent to diffusive 
propagation occurs when (V/AL) =V/ (a r/2) ~ 1, or when r -2 (V/A), we 
assumed that the injection of particles occurred at z=2(V/A), rather than at 
z = 0, and at time t = z/V + = ( ~ ~ ), which is the time it takes the 

supercoherent pulse to reach the transition point at speed V . 

The function (z,t) was obtained by including only the first two 
terms in the eigenfunction series that represents the particle distribution 
function. If the next two higher order terms are included, the 6- function 
pulses in the density assume a finite spread due to dispersion. Physically, 
dispersion describes statistical fluctuations within a bunch of particles 
whose velocities wander randomly in direction over one hemisphere with an 
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average speed along the field of approximately 0.6V. An approximate 
solution of this more rigorous treatment, valid when (V/AL) £ 1, is 
obtained by convolving F 0 + { z ,t} with a coherent Gaussian (Earl 1976b), 


F 0 O*,t’> = 


r 


exp{- (z-z n ) / 4D^t] 
dz -2(vD *t)% “ 


F 0 { V t} 9 


(5) 


2 

where D* = 0.0405 (V /A) is the coefficient of dispersion introduced by 
Earl (1974a, eq. 83) and discussed by Kunstmann and Alpers (1977), and where 
z' = z + 2 (V/A), f = t + ( f + }( l~). 

Figure 1 shows as solid lines the density profiles for 40 MeV protons 
predicted by equation (5) for (V/AL)=0.2 and 0.7. To illustrate the 
importance of the coherent effect, the dashed lines .show separately the 
Gaussian profiles into which the coherent 6 -functions are converted by 
the convolution operation. At the right, where the choice (V/AL) =0.7 
leads to virtually coherent or scatter-free propagation, the coherent 
Gaussian is a dominant feature which overshadows the inconspicous tail 
that arises from the diffusive wake. 'But, even at the left, where the 
choice (V/AL) =0.2 leads to relatively diffusive propagation dominated 1 
by the wake, the peak intensity provided by the coherent Gaussian amounts 

to a substantial fraction (36%) of the total intensity present, 
at' the time of maximum. To test the accuracy of ; the 

approximations leading to -equation (5), John.Bieber (1977) calculated 
density profiles which apply to exactly the same conditions assumed in 
Figure 1, but tcek rigorously into account' all of the first four 
eigenfunctions . Evidently, the differences between the rigorous calcul- 
ation, whose predictions appear as dotted lines in Figure 1, and the 
approximate treatment invoked here are insignificant. 
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Equation (5) gives the particle density following an impulsive 
injection at the sun. On comparing the temporal profiles of de'nsity 
calculated at z’ = 1.2 AU with the temporal profiles of western flare 
events, we found that, in almost all cases, the calculated density 
successfully fits the events during the decay phase, but the initial 
rise to maximum is consistently shorter than that observed. Although 
perpendicular diffusion might cause this discrepancy, this possibility 
is considered unlikely because it would imply that the mean propagation 
time x has a form (Jokipii and Parker 1968, 1969) 


T “ K “ 6 
that is not obseved. 

A second, more viable, possibility is that the additional delay "is 

due to a coronal effect which makes the release of. particles from the sun 
persist over a finite time interval. (Feit 1973; Nielsen et al. 

1975; Palmer et al. 1975; Reinhard 1975; Ng and Gleeson 1976; 

Wibbexenz et al. 1976). Because equation-' (5) refers to an 
injection whose temporal profile can be represented by a 6 -function,, the 
solution for an injection f(t) of finite duration is the convolution of 
f(t) with this Green's- function. As. an approximation, it was assumed that 
the injection profile i's a Gaussian with a standard deviation a. The- peak 
of the Gaussian was placed at t = 3a and the Gaussian was truncated: in such 
a way that particles are not released before t = £F. The convolution of the 


Gaussian with equation (5) is: 
F{z',t} = f 


0 


dt- F 0 {z’,f} 




( t-t '-3 j) ' 
■ 2<j' a 


C2tt)\ 


( 6 ) 


where t is the time of observation. The lower limit of this integral is 
zero, because F 0 {z',t f } = 0 for f < 0. The upper limit is t, because the 
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Gaussian was truncated for t-t T < 0. Clearly, this is not necessarily the 
only possible solution ^tfhich can describe a particle event. However, it will 
become evident that, for most events, this injection profile does give a 
satisfactory representation of the initial phase. 
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IV. MODEL FITTING 


The entire time-intensity profiles of thirty events were fitted 
by- -eye -using t-hi-s model . To demonstrate the fitting of a typical 
event, such as the one that occurred on July 7, 1970, we show in 
figure 2 the 0. 5-1.1 MeV electron and 30-50 MeV proton time-intensity 
profiles compared with densities calculated from equation (6). Solid 
curves give profiles for the values of (V/AL) which best fit the observed 
particle density during the decay phase. Dashed curves give profiles 
for two neighboring values of (V/AL). After the most appropriate value 
of (V/AL) was found by comparing the calculated and observed decay 
phases, the Gaussian injection width o was varied in steps of 0.2 hours 
until a best fit over the entire profile was found. For the event in 
Figure 2, this procedure leads to the following estimates of (V/AL) and 
o for electron and n. 40 MeV protons j 
(V/AL) = 0.06 ±0.005, 


a = 0 . 4 ±0 . 1 hours 
(V/AL) = 0.14 ±0.01, 


0. 5-1.1 MeV e 


30-50 MeV p 

a = 0.6 ±0.1 hours. 

In figure 2, note that the rigidity of the electrons is 2^ orders 
of magnitude smaller than that of the protons, but that the scattering length 
of the electrons is only a factor or 2 smaller than that of 40 MeV protons. 
This weakness of rigidity dependence in interplanetary propagation has been 
found in many events. A similar conclusion was drawn by Countee and 
Lanzerotti (1976). As will be discussed later, the rigidity dependence of 
the coronal injection profile is even weaker. 


Figure 3 shows a sample of 0.5-1 .1 MeV electron events. Figure 4 shows 
a corresponding sample of 40 MeV proton events. These events cover the 
entire range of (V/AL) from coherent propagation to diffusion. In each case. 
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to allow proper study o£ the behavior near onset, the background flux was 
subtracted. For purposes of comparison, the density calculated for the 
indicated parameters was averaged over the same time interval as the data. 
Figures 3 a and 4 a show the time histories for an event on December 29, 

1967 (W77°) which had the extremely sharp time profile of scatter-free 
microevents. Because this event was very brief, electron data were sampled 
every three minutes, and proton data were averaged over fifteen-minute 
intervals. -In this case, the best fit was obtained fox o = 0 which corresponds 
to impulsive injection. The ratios (V/AL) = 0.48 for 0. 5-1.1 MeV electrons, 
and (V/AL) = 0,85 for 40 MeV protons obtained for this event are the highest 
of all events -analyzed. These large values of (V/AL) correspond to highly 
coherent propagation. However, note that the propagation- of protons was . 
more coherent than that of electrons. Figures 3 b and 4 b show a microevent' 
on August 3, 1967 (W85 ) during which electron and proton data were averaged 

oyer fifteen minutes and thirty minutes, respectively. Essentially the same • 
injection width, o -=;0,4 - hr , describes the initial phase of both electrons 
and 40 MeV protons. 

The remaining events in Figures 3 and 4 become progressively more 
diffusive as (V/AL) decreases. These include the 0.5-1. 1 MeV electron and 
a- 40 MeV proton events on August 14, 1968 (W80°); February 24, 1969 (W32°) ; 
February 25, 196? (W37°) ; April 6, 1971 (W80°) November 24, 1969 (W31°) ; 
and March 27, 1969 (W68°). Figure 3 g shows, in addition, the 3-12 MeV 
electron .event on February 25, 1969 (W37°). This event was among the very 
few for which 3-12 .MeV electron flux could be analyzed. In general, the 
electron flux at 3-42 MeV was masked by a combination of a high electron 
background and a steep electron spectrum. By comparing the electron time 
histories in the interval of 0.5-1. 1 MeV .and 3-12 MeV (Figures 3 d and 3 g) , 
we note xhat electrons have essentially the same injection profile and 
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propagate in about the same manner over the velocity range 6 = 0. 9-1.0 
and rigidity range R = 1-8 MV. The corresponding 40 MeV protons (the 

second event in Figure 4 d) are characterized by a value of (V/AL) which 

is- only a factor of 2,5 larger. 

For the event of April 6, 1971 (W80°) , which is a typical western 
hemisphere event. Figure 5 compares the propagation characteristics of 
protons over the range from 4 to 80 MeV, corresponding to a range of 

rigidities R from 87 to 400 MV, and a range of velocities g from 0.09 to 

0,39. Figure 3 e shows the temporal profile of the 0. 5-1.1 MeV electrons 
for this event. Solid curves give the best estimates of (V/AL) in each 
range, while dashed curves show possible fits in steps of 0.02 about the 
best estimates of (V/AL) , Here (V/AL) increases with decreasing proton 
energy. This means that propagation of lower energy protons is more coherent 
than that .of higher energy ones. Since the mean free path is proportional 
to (V/AL), for a fixed value of L, this result also means that lower energy 
protons have a larger mean free path. On the other hand, comparison of 
Figures 3 e and 5 shows that electrons propagate more diffusively than 
protons. Figure 5 also shows that the injection profile broadens as the 
velocity of the particles decreases. 

Although the model gave a good account of the time profiles for 30 
out of a total of 31 western hemisphere events, a single exception occurred 
on March 12, 1969. The identification of the flare was based on strong 
associations with Ha, X-ray emission, type II, III, and IV radio bursts, and 
intense microwave bursts. But the proton profiles between 20 and 80 MeV, 
shown in Figure 6, deviate significantly from predicted profiles given by 
the dashed curves. The best estimates of (V/AL) indicate that protons at 
higher energies are more coherent than those at lower energies. This feature, 
which disagrees with the behavior of other events, has no satisfactory 
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explanation at the present time. 

The analysis shows that the theory of focused diffusion combined 
with a Gaussian injection profile accounts satisfactorily for the time 
histories of both proton and electron events. The ratio (V/AL) is, on 
the one hand, smaller for electrons than for protons, and on the other 
hand, larger for lower energy protons than for higher energy ones. This 
result has an important bearing on the propagation of low energy particles 
in interplanetary space. The model is especially successful in describing 
the scatter- free micro events, characterized by a very sharp time profile 
with a rapid rise and decay (Lin et al. 1973; Lin 1974; McDonald and Van 
Hollebeke 1973), as shown in Figures 3 and 4. Classical diffusion theory 
is totally inadequate in explaining these events-, but they axe success- 
fully explained by the theory of focused diffusion, indicating that the 
coherent effect dominates over scattering in such cases. 
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V. RESULTS AND DISCUSSION 


Detailed information on (V/AL) and the Gaussian injection width 
for 3CT events from ^W20° toMtflOO 0 is given in Table I. Figure 7 shows, 
for 18 events, the frequency distribution of (V/AL) for 0. 5-1.1 MeV 
electrons, 30-50 MeV protons, and 9-19 MeV protons. Excluded from this 
distribution are events for which the analysis could not be performed 
simultaneously for all three of these particle categories. It is 
apparent that (V/AL) is generally smaller for electrons than it is for 
protons. Similarly, for 9-19 MeV protons, the distribution, which extends 
to relatively large values of (V/AL), is broader than it is for 30-50 MeV 
protons. 

Figure 8 shows six examples of the mean free- path A, calculated from 
equation (4) for individual events, as a function of rigidity for electrons 
and protons. Here, data at R - 100 MV refer to protons, while those at 
R =1 and 8MV refer to 0.5-1. I and 3-12 MeV electrons, respectively. In 
two case&, the events on September 29-, 1968 and November 24, 1969, proton 
data were extended to 326 MeV (rigidities up. to 847 MV) using an essentially 
identical detector on. 0G0 V which was described by Jones et al. (1967). 
Clearly, the rigidity dependence of A fox protons varies from event to event, 
but, between 100 and. 800 MV, the mean free path either stays relatively 
constant or increases with decreasing proton rigidity. This behavior is 
also exhibited in the average A over all events shown in Figure 9, which 
includes all events except the microevents. The micxoevents were excluded 
because their limited number and incomplete analysis would create a bias. 

The error bars in Figure 8 give the standard deviation of the average. 

Note that for protons between 4 and 80 MeV (R = 100-400 MV) , A increases 
slowly with decreasing rigidity. For protons in this range, a typical mean 
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free path is ^0.2 AU. The 1 MV electrons have a mean free path approxi- 
mately 2 to 3 times smaller. 

The inverse relationship between the mean free path and particle rigidity, 
idiich has also been reported by Wibberenz [1977) and Zwickl 

(1976) , is in strong- disagreement with quasilinear theory of pitch-angle scatter- 
ing. It follows from equations (2) and (4) above that A « ^ (Jokipii 1971) . 

Thus, for q < 2, the quasilinear theory predicts, in disagreement with the 
observed dependence, that the mean free path decreases with decreasing rigidity. 
On the other hand, if the experimental result that X either stays constant 
or increases with decreasing proton rigidity is interpreted within the context 
of quasilinear theory to mean that q >_ 2, then purely coherent propagation 
(Earl 1974a) is predicted but not observed. Moreover, the diffusion coeffi- 
cient calculated fallowing Jokipii and Coleman (1968) from the measured in- 
terplanetary magnetic field spectrum is approximately one order of magnitude 

22 2 

smaller than the observed value of ^(0.6-1.2)xl0 cm /sec at 40 MeV. 

The paradox is further compounded at low rigidities by the electron 
data, which do not conform to the monotonic inverse dependence found fox 

protons, but which suggest instead that the scattering length is a compli- 
cated function of rigidity which goes through a maximum somewhere between 
10 and 500 MV, and then dips to minimum somewhere above 1000 MV. This 
implication must be taken very seriously, for it follows from. the nature 
of the Eorentz force that the scattering length, idiich characterizes the 
effect of a static magnetic field, must be a function of rigidity alone. 

In the absence of a definitive theory, a complicated dependence of 
pitch-angle scattering upon rigidity cannot be ruled out a priori. Never- 
theless, because this complexity seems implausible, it is worth considering 
alternatives that retain a simple dependence. One possibility is that both 
of the parameters, q and A, which specify <f> in equation (2), are functions 
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of rigidity. This hypothesis, that the anisotropy of scattering, as well 
as the amplitude, depends upon rigidity, means that the scattering length 
(V/A) for a given mean free path A increases as q decreases (see eq. 4). 
Figures 8 and 9 show that the mean free path of electrons is smaller 
than that of protons. If q is held fixed at q = 1.5, as it was in’ the 
above analysis, this also means that (V/A) is smaller fox electrons than 
protons. But if scattering at the low rigidities of the electrons is more 
isotropic than it is at the rigidities of the protons, which are two 
orders of magnitude larger, then two different mean free paths could 

i 

correspond to the same scattering length. For example, if q = 1.0 for 
electrons. (i,e. A =V/A) and if q = 1.5 for protons (i.e. A = 2.4 V/A) , 
then the electron mean free path, for a given (V/A), is a factor of 
2,4 smaller than that of protons. This number is close enough to the 
observed ratio to make tenable an interpretation of the maximum in the A 
vs. rigidity curve as the effect of a slow increase of q superposed upon 
a slow decrease of (V/A) . However, this interpretation is not very satis- 
fying, for it offers no fundamental explanation of the postulated rigidity 
dependences. 

To seek further alternatives that embody monotonic dependences, it is 
appropriate to express the results appearing in Figure 9 as the following 
power law in the two particle parameters intrinsic to electromagnetic 
forces, the velocity V and the charge to mass ratio (Q/m) : 


V 

A 


C S ) 

v m ' 


- 0.1 


-0.3 

V 


In this expression, the dependence on (Q/m) is weak enough to suggest a 
scattering process that depends on velocity alone. If the relationship 
R = (mV/Q) is invoked to convert to a power law in (Q/m) and rigidity R, 
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the result is : 


V 

A 


ec 



3 


CQ/m) 


-0.4 


Because scattering by static magnetic fields must depend on rigidity alone, 
the presence here of a factor in CQ/ m ) demonstrates again the conclusion 
reached earlier that scattering is not a simple function of rigidity. 
Although the presence of electrical forces in interplanetary space seems 
very unlikely, they constitute the only other possible influence on charged 
particle trajectories. To investigate this possibility, we express the 
power law in terms of the kinetic energy per unit charge (E/Q) , 


j -CE/Q) 


-0.15 


CQ/m) 


-0.25 


Here again, the dependence embodies a factor in (Q/m) which is not consistent 
with the pure (E/Q) dependence expected if the scattering is produced by 
electric fields alone. These considerations lead back to the conclusion 
that a non-monotonic dependence of scattering length on rigidity is the most 
plausible explanation of our results. 

Many investigators have analyzed propagation in terms of the distance 
travelled by particles from the time of flare acceleration. For example, 
Bryant et al. (1965) found that, during the event on September 28, 1961, 2-300 

MeV protons shared the same probability distribution in distance travelled, 

and concluded that propagation during this event was rigidity independent 
(i.e. K^BR 0 ). Reinbard and Wibberenz (1974), and Ma Sung et al. (1975) 

showed that, at the time of maximum intensity, distance travelled could be 

expressed as a linear function of velocity. 


VAT = C + B V 
max max max ’ 


(7) 


where A T ^ = risetime to maximum intensity, where C max was interpreted 
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as the distance travelled in interplanetary space and where B max was 
interpreted as the time spent in the corona. This interpretation is 
equivalent to the conclusion that the diffusion coefficient K is independent 
of rigidity. The present analysis indicates that the dependence of propagation 
upon rigidity. is weak, but not entirely negligible. Furthermore, it is 
not possible to separate the injection and interplanetary effects^ as 
was done in equation (7), because the coronal injection profile also displays 
a weak dependence on velocity. 

In section IV, it was shown that the injection width o increases with 
decreasing particle velocity. For example, in Figure 3, for the event 
on April 6, 1971, cr =0.2 hour for electrons (3 = 0.91). For protons, in 
Figure 5, a =0.3 hours at SO-SO MeV (3 = 0.35), and increases to 3.0 hour 
at 6-9 'MeV ( 3 =0.125). This is also evident in the frequency distribution 
of a shown in Figure 10, which is based on the same set of 18 events as 
Figure 7, The majority of events have an injection width o hour, but a 

few are characterized by a large a of several hours. For those events with 
>1 hour at 40 MeV, it was found that a Gaussian profile truncated at the 
3cr level gives a slower onset phase, than is observed. A Gaussian truncated 
at the 1.50 level appears to be more appropriate, but this change does not 
affect events with a narrow injection profile. Figure 11 shows five proton 
events with prolonged injections each in a different energy range. Injection 
widths up to 8 hours have been found for proton events. The finite duration 
of injection can only arise in the corona. However, the methods employed here 
do not distinguish between prolonged acceleration and continuous injection. 
Since no apparent correlation is found between the injection width and the 
duration of the type IV bursts, it is assumed hereafter that a characterizes 
the duration of particle release from the corona. 

Figure 12 shows a plot of a versus particle velocity averaged over- all 
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events with a $1 hour at 40 MeV. Events with prolonged injection were 
excluded because their analysis was incomplete. The large magnitude of 
the error bars, which indicate the standard deviation of the average, 
arose from large variation from event to event, and not from the un- 
certainties in individual measurements of a. Data at $ > 0.9 refer to 
electrons, and the rest refer to protons. There appears to be a continuous 
variation of o as a function of the velocity, regardless of the rigidity. 

Because knowledge of the coronal release process is very limited, we 
assume, for simplicity, that the injection width follows a power law in 
rigidity and velocity. Then the best fit to the data leads to 

cr « r- 0. 07±0. 06 g-0.55±0.85 

in which the rigidity dependence is negligible. The dashed curve in 
Figure 12 shows a fit for the data given by the relation: a 0-55. 

Figure 13 shows two typical events, on September 29, 1968, and February 
25, 1969, for which the relation given above appears to be appropriate. 
However, deviations do occur. Figure 14 shows two examples. During the 
event on April 22, 1971, the injection width was constant over the range 
0.1 t B - 0.9. During the event on March 29, 1970, a exhibited a strong 
dependence on £5, for a = 0 for 3-12 MeV electrons (3 = 1.0), and a increases 
six-fold as proton velocity decreases by a factor of three. However^ this 
eyent was not included in the average shown in Figure 12, because a > 1 
hour at 40 MeV. These examples illustrate variability in the coronal 
release. However, one feature of the coronal release remains unchanged 
throughout all events: the electrons have an injection width similar to that 
of protons above 40 MeV, even though their rigidity is 2h orders of magnitude 
lower; on the other hand, low energy protons are generally associated with 
broader injection profiles than high energy protons. These tendencies 
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suggest strongly that rigidity plays essentially no role in the release 
process. 

The velocity dependent injection width implies that the mean escape 
rated also depends on velocity: 


2. 35 o 


0.55 

3 


where we approximated the duration of' injection by the full width at half 
maximum of the Gaussian. This relation predicts that high energy proton 
fluxes in space should have a steeper longitudinal gradient than low energy 
ones. The absence of rigidity as a factor in coronal release severely limits the 
possible mechanisms for this process, as discussed by Newkirk and Wentzel (1978) . 
In a separate paper by Ma Sung, Van Hollebeke, and McDonald (1978), this 
question will be examined by analyzing the characteristics of eastern 
hemisphere solar events. 

Because both coronal escape and interplanetary propagation depend on 
particle parameters, the energy spectrum measured at the time of maximum 
intensity at a given position in space does not necessarily reproduce the 
source spectrum. Interplanetary propagation tends to enhance fluxes at 
lower energies, because propagation of low energy particles is more coherent. 

On the other hand, the velocity dependent injection profile tends to 
enhance fluxes at higher energies. Neither effect is strongly dependent 
on particle parameters, and the two effects work in opposite directions. 

On the basis of a careful study of many events, we conclude that the com- 
bined process alters the source spectrum by no more than 5% at 1 AU. 

This analysis provides evidence that the average proton energy spectrum 
that was measured by Van Hollebeke et al. (1975), for events in W20° to 
W80°, is identical to the source spectrum to within 5%. 
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Because the predictions of the model accurately describe the entire 
evolution of particle events from the western hemisphere, the time of 
particle release from the sun can be determined from the model, as listed 
in column 7 of Table I, By comparing this release time with the onset of 
the 2-20 X soft X-ray bursts (listed in column 3 of Table I), we found that, 
with only one exception,- no. delay occurred between the predicted particle 
release time and the soft X-ray onset, provided that the onset time of 
X-ray emission closely coincides with the onset time of type II and type IV 
radio bursts or intense microwave bursts. In the single exception, on June 
29, 1971, the predicted release of protons occurred 3 hours after the 
occurrence of X-ray emission and type II and IV bursts. Four other events 
show significant delay in particle release compared to X-ray onset. Of 
these, three events, on January 28, 1970, January 29, 1970, and April 22, 1971, 
were associated with neither spectral radio bursts nor microwave bursts. 

Because flare assocation for these events was weak, the delays are probably 
not meaningful. The fourth event, on March 23, 1970, shows a delay of 2 
hours if the predicted particle release time is compared with the X-ray 
onset, but shows no delay if it is. compared with the onset of type IV hursts, 
which occurred 2 hours after the onset of soft X-ray emission. Thus, in 
almost all cases, particle release begins at flare acceleration indicated by 
the concurrent outbursts of Ha emission, soft X-ray emission, type II and IV 
bursts, and microwave bursts. 



26 


VI. SUMMARY AND CONCLUSION 

In this paper, a propagation model, which combines a Gaussian 
profile for particle release from the sun, with interplanetary particle 
densities predicted by the theory of focused diffusion, meets with 
remarkable success in describing the time -intensity profiles of 30 flare 
associated proton and electron events from the western hemisphere of the 
sun. This model makes possible a clear separation_of coronal and inter- 
planetary effects. By comparing the predictions of the model with observed 
events over a wide range of propagation conditions and particle rigidities, 
we obtained the results- summari-zed in the .next' two -paragraphs , 

The scattering mean free path X for flare associated particLes is large. 
Its magnitude is on the order of 0.1 to 0.3 AU for 4-80 MeV protons, and 
a factor of 2 to 3 smaller for 0. 5-1.1 MeV and 3-12 MeV electrons. Because 
the rigidities here differ by more than two orders of magnitude, this result 
indicates that \ is only weakly dependent on rigidity. The mean free path 
is at least one order of magnitude larger than is predicted by the quasilinear 
-theory of pitch-angle scattering. Furthermore, X exhibits a complex 
dependence on rigidity that is in strong disagreement with quasilinear 
theory. 

The release of particles into interplanetary space’ generally occurs 
over a finite time interval. If the injection profile is approximated by 
a truncated Gaussian, then the r.m.s, width cr is less than 1 hour except 
for a few events in which it can be as large as several hours. The injec- 
tion width decreases with velocity, but it does not depend on rigidity. In 
almost all events, the release of particles occurs at the time of flare 
acceleration Indicated by concurrent events in Ha emission. X-ray emission. 



27 


type II, III, and IV radio bursts, and intense microwave bursts. 
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Table I. SOLAR PARTICLE EVENTS - INTERPLANETARY PROPAGATION 

MODEL FITTING 




ACCEL. 




T 




TIME 

V 

V 


T inj 

PARTICLE 

DATE 

LONGITUDE 

CUT) 

.(AU/HR) 

AL 

a (HR) 

(UT) 

SPECIES 

5/28/67 

1732 

5.35 

6.55 

- 



e 




2.52 

0.10 

0 

5. 9-0. 5 

P 




1.98 

0. 12 

0.15 

5.4±0.S 

P 




1.58 

0.12 

0 

5.8±0.5 

P 




1.15 

0.14 

0.4 

5.240.5 

P 




0.79 

0.16 

0.8 

4.5540.5 

P 

8/3/67 

W85 

9.13 

6.55 

0.14 

0.3 

9.240.25 

e 




1.98 

0.46 

0.4 

9.5±0.5 

P 




1.58 

0,50 

0.4 

9. 3+0. 5 

P 

12/3/67 

W105 

8.78 

6.55 

0.03 

0.2 

8.2+0.25 

e 




2.52 

0.08 

0.4 

8.9±0.5 

P 




1.98 

0.10 

0.4 

9. 4+0. 5 

P 




1.58 

0.10 

0.4 

8.9±0.5 

P 




1.15 

0.12 

0.6 

8.3±0.5 

P 




0.79 

0.16 

0.6 

8.5540.5 

P 

12/29/67 

1V77 

0.63 

6.55 

0.48 

0 

0.9540.1 

e 




1.98 

0.85 

0 

0.9540.25 

P 


_ 


1.58 

- 



P 




1.15 

a.0.75 

0 

1.4±0.25 

P 




0.79 

^0.75 

0 

2.1±0.25 

P 

2/17/67 

W47 

2.72 

6.55 

- 



e 




2.52. 

Q. 18 

0-4 

2.340.. 5 

• P 




1..98 

0.18 

0.4 

2.3±0.5 

P 




1.58 

0.18 

0.4 

2.2540.5 

P 




1.15 

0.20 

0.6 

2.1±0.5 

P 




0.79 

0.24- 

0.6 

2. 3+0. 5 

P 





0.26 




8/14/68 

W80 

13.28 

6.55 

0.09 

0.6 

12.7 40.25 

e 




2.52 

0.20 

0.8- 

13.040.75 

P 






1.0 






1 -98 

0.20 

0.8- 

12.840.5 

P 






1.0 






1.58 

0.20 

0.8 

12.7±0.5 

V 




1.15 

0.30 

1.0 

12.4±0.5 

P 




0.79 

- 



P 


The error in (V/AL) is ±0.01, and the error in a is ±0.1 hr. 

o 

The acceleration time is assumed to he at the onset of 2-20 A X-ray burst. 
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TABLE I. (Con't) 




ACCEL. 




T. . 



TIME 

V 

V 


in] 

DATE 

LONGITUDE 

CUT) 

(AU/HR) 

AL 

cf(HR) 

(UT) 

9/29/68 

W5 2 

16.15 

6.55 

0.05 

0.4 

15. 7± 0.25 




4.68 

0.10- 

0.12 

0.4 

15.9±1 .0 




4.32 

0.10- 

0.12 

0.4 

15.9±1 .0 




3.74 

0.10- 

0.12 

0.6 

15 . 9± 1 . 0 




3.17 

0.12 

0.6 

15.9+1.0 




2.52 

0.12- 

0.14 

0.6 

15 . 9±0 .5 




1.98 

0.12- 

0.14 

0.8 

15.9+0.5 




1.58 

0.12- 

0.14 

0.8 

,15.8±0.5 




1.15 

- 






0.79 

- 



10/4/68 

W37 

(03)23.62 

7.20 

0.06 

0.6 

24 . 0± 0 . 5 




6.55 

0.06 

0.6 

23.6+0.25 




2.52 

0.12 

0.6 

23.9±0.5 




1.98 

0.12 

0.6 

23.9+0.5 




1.58 

0.12 

0.8 

23.8+0.5 




1.15 

0.12 

1.0 

23.3±0.5 




0.79 

- 



10/31/68 

W37 

(30)23.52 

1.58 

0.20 

^6 . 0 

1. 0+1.0 

2/24/69 

W32 

22.9 

6.55 

0.08 

o;3 

22.5±0.4 




1.98 

0.20 

0.4 

22.3±0.5 




1.58 

0.20 

0.6 

22.7+0.5 




1.15 

0.20- 

0.22 

0.6 

22,1+0.5 




0.79 

0.24 

0.6 

22.8±0.5 

2/25/69 

W37 

8.80 

7.20 

0.07 

0.3 

8.7±0.5 




6.55 

0.07 

0.3 

8 . 7±0 .4 




2.52 

0.16- 

0.18 

0.4 

8. 8+0. 5 




1.98 

0.16- 

0.18 

0.6 

8. 5+0. 5 




1.58 

0.16- 

6.18 

0.8 

8.25+0.5 




1.15 

0.18- 

0.20 

0.8 

8. 1±0 .5 


PARTICLE 

SPECIES 

e 

P 

P 

P 


P 

e 

e 

P 

P 

P 

P 


0.79 


P 


T3’-d^3tt> >d t) T! d d d (P <( d >d d d ■d’d 
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TABLE I. (Con't) 


DATE 

LONGITUDE 

ACCEL. 

TIME 

CUT] 

V 

(AU/HR] 

2/27/69 

W65 

13.82" 

6.55 

2.52 

1.98 

1.58 

1.15 

0.79 

3/27/69 

W68 

13.18 

6.55 

2.52 

1.98 

1.58 

1.15 

0.79 

9/25/69 

W14 

6.50 

6.55 

2.52 

1.98 

1.58 

1.15 

0.90 

11/24/69 

W31 

9.17 

6.55 

3.74 

3.17 

2.52 

1.98 




1.58 

1.15 




0.90 

12/30/69 

W85 

18.92 

2.52 

1.98 

1.58 

1/28/70 

W26 

10.22 

6.55 

2.52 

1.98 

1.58 

1.15 

0.90 




T. . 


V 


mj 

PARTICLE 

AL 

a (HR] 

CUT] 

SPECIES 

0.02 

0.8 

13.5±0.4 

e 

0.12 

2.0 

13.9+0.5 

P 

0.16 

3.0 

13.840.5 

P 

0.20 

3.0 

14.5+0.5 

P 

0.26 

4.0 

14.5±0.5 

P 

- 



P 

0.02 

0.2 

13.140.25 

e 

0.06 

0.8 

14,4+0.5 

P 

0.06 

0.8 

13.940.5 

P 

0.06 

0.8 

13.9+0.5 

P 

0.06 

0.6 

13.940.5 

P 

- 



P 

0.07 

0.2 

7.040.25 

e 

0.20 

0.4 

6.840.5 

P 

0,20 

0.6 

6.8+0.5 

P 

0.20 

0.6 

6.740.5 

P 

0.22 

0.6 

6.140.5 

P 

0.22 

0.6 

6.540.5 

P 

0.02 

0.2 

8.040.25 

e 

0.06 

*0.2 

9.0+1.0 

P 

0.06 

*0.2 

9.041.0 

P 

0.06 

0.2 

8.940.5 

P 

0.06- 

0.2 

8.940.5 

P 

0.08 




0.08 

0.2 

8.940.5 

P 

0.10- 

0.6 

9.840.5 

P 

0.12 




0. 12- 

0.6 

9.740.5 

P 

0.14 




0.12 

2.0 

19.4+0.5 

P 

0.12 

3.0 

18.9+0.5 

P 

0.12 

2.0 

19.540.5 

P 

0.04 

0.2 

12.540.25 

e 

0.18 

1.0 

13.8+0.5 

P 

- 



P 

0.22 

2.0 

14.2+0.5 

P 

- 


- 

P 

- 



P 
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TABLE I . (Cont ' d) 




ACCEL. 




T. . 




TINE 

V 

V 


inj 

PARTICLE 

DATE 

LONGITUDE 

CUT) 

(AU/HR) 

AL 

a (HR) 

(UT) 

SPECIES 

1/29/70 

W42 

10.17 

6.55 

- 



e 




2.52 

0.07 

2.0 

11.9+0.5 

P 




1.98 

0.12 

4.0 

11.9+1.0 

P 




1.58 

0.12 

4.0 

12.5*1.0 

P 




1.15 

0.16 

4.0- 

13.6+1.5 

P 






5.0 






0.90 

0.18 

- 


P 

1/31/70 

W62 

14.98 

6.55 

0.03 

0.6 

14.2±0.25 

e 



2.52 

0.10 

0.6 

14.4*0.5 

P 




1.98 

0.10 

0.6 

14.4*0.5 

P 




1.58 

0.12 

0.8 

14.3*0.5 

P 




1.15 

0.14 

1.2 

14.2+0.5 

P 




0.90 

- 



P 

3/23/70 

W62 

15.58 

6.55 

0.05 

0.4 

17.2*0.25 

e 



2.52 

0.10 

0.2 

17.9+0.5 

P 




1,98 

0.12 

0.2 

18.1*0.5 

P 




1.58 

0.12 

0.2 

17.8*0.5 

P 




1.15 

0.14 

0.4 

17.2*0.5 

P 




0.90 

0.16 

0.4 

17.6+0.5 

P 

3/29/70 

W37 

0.37 

7.20 

0.015 

0 

1.0+0. 5 

e 




6.55 

0.015 

- 


e 




2.52 

0.04 

1.0 

0.0*0. 5 

P 




1.98 

0.05 

2.0 

0. 0*0.5 

P 




1.58 

0.06 

2.0 

0.4*0. 5 

P 




1.15 

0.10 

5.0 

0.2*0. 5 

P 




0.90 

0.12 

6.0 

0.2+0. 5 

P 

7/6/70 

W90 

21.45 

6.55 

0.06 

0.3 

21.0*0.25 

e 



2.52 

0.14 

0.2 

21.9*0.5 

P 




1.98 

0.16 

0.6 

21.3*0.5 

P 




1.58 

0.16 

0.8 

21.3*0.5 

P 




1.15 

0.16 

0.6 

21.2*0.5 

P 




0.90 

0.18 

0.8 

21.6*0.5 

P 

7/7/70 

W100 

16.70 

6.55 

0.06 

0.4 

16.5+0.25 

e 



2.52 

0.12 

0.4 

16.9*0.5 

P 




1.98 

0.14 

0.6 

16.8*0.5 

P 




1.58 

0.16 

0.8 

16.8*0.5 

P 




1.15 

- 



P 




0.90 

- 



P 

1/24/71 

W49 

23.0 

1.15 

0.18 

8.0 

23.5*0.5 

P 
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TABLE I. CCont’d) 




ACCEL . 




T 



TIME 

V 

V 


inj 

DATE 

LONGITUDE 

CUT) 

CAU/HR) 

AL 

a (HR) 

CUT) 

4/6/71 

W80 

9.42 

6. 55 

0.04 

0.2 

9.4±0.4 




2.52 

0.12 

0.3 

9.910.5 




1.98 

0.14 

0.4 

9.810.5 




1.58 

0.16- 

1.4 

9.210.5 





0.18 






1.15 

0.18 

1.0 

9.910.5 




0.90 

0.20 

> 1.0 

10.010.5 

4/20/71 

W50 

19.32 

6.55 

0.06 

0.3 

19.010.25 




2.52 

- 






1.98 

0.20 

0.6 

19.310.5 




1.58 

- 






1.15 

0.22 

0.6 

20.110.5 




0.9 

- 



4/22/71 

W61 

%n 

CO 
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Figure 1 


Figure 2 


Figure 3 


FIGURE CAPTIONS 


Comparison of density profiles for 40 MeV protons (V = 1.98 AU/hr)i . 

The solid curves show the densities according to* the convolution 
approximation given by equation (5) , the dotted curves show those 
according to the accurate solution calculated by Bieber (1977), and 
the dashed curves show the contribution of t-he coherent Gaussian. The 
difference between the solid and dotted curves is not significant 
for both (V/AL) = 0.2 and 0.7. The contribution of the coherent 
Gaussian overshadows the inconspicuous tail of the diffusive wake 
when (V/AL) = 0.7, and it amounts to a significant fraction (36%), of 
the total density at maximum even when (V/AL) = 0.2 for which the 
diffusive wake dominates. 

Comparison of the observed time intensity profiles of 0. 5-1.1 MeV 
electrons and 30-50 MeV protons with the calculated density profiles. 
The solid' curves give profiles for the- values of (V/AL) which best 
fit the observed particle density during the decay phases. The dashed- 
curves give profiles for two neighboring values of (V/AL) . The value 
of a indicates the standard deviation of the Gaussian coronal injection, 
profile which best fits the onset -to -maximum phase of the event. 

A sample of the 0. 5-1.1 MeV electron events fitted by the calculated 
density profiles. The electron data were sampled every 3 minutes in 
(a), and were averaged over 15 minutes in the rest. The sample shows: 
(a) and (b) , two scatter-free electron events with large (V/AL); (c). 
to (£) , progressively more diffusive events with decreasing (V/AL)-; 

(g), an electron event at 3.-12 MeV. The same event at 0. 5-1.1 MeV is 
shorn as the second; event in (d) . 



37 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9 


A sample of the 40 MeV proton events fitted by calculated 
densities. The proton data were averaged over 15 minutes in (a) , 
and over 30 minutes in the rest . The- sample. shows , (a) and (b) , 
two scat-tex-free proton events with large (V/AL ) ; (c) to (g) , 
progressively more diffusive events with decreasing (V/AL) .. 

Comparison of the propagation characteristics of protons between 
4 and 80- MeV (R = 87-400 MV) for a typical western event. The solid 
curves give the best estimate of (V/AL) in each energy range, and 
the dashed curves show the variations about the best fit in steps of 
0.02 in (V/AL). Note that (V/AL) increases with decreasing proton 
energy. 

Proton time intensity profiles for the March 12, 1969 event. The 
dashed curves show the best possible fit according to Equation (6) . The 
fit is poor for > 20 MeV protons. There is no satisfactory explanation 
for this behavior. 

Frequency distribution of (V/AL) for 18 0. 5-1.1 MeV electron, 30-50 
MeV proton, and 9-19 MeV proton events. For electrons, (V/AL) is generally 
much smaller than it is for protons. The lower energy protons tend 
to have a broader distribution in (V/AL) than the higher energy ones do. 

The scattering mean free path A. as a function of rigidity R for six 
electron and proton events. Data at R - 100 MV refer to protons above 
4 MeV, and those at 1 and 8 MV refer to 0. 5-1.1' and 3-12 MeV electrons, 
respectively. • Between 100 and 800 MV the proton mean free path either 
stays constant or decreases with rigidity. The electron mean free path 
is generally a factor of 2 to 3 smaller than that of protons above 20 MeV. 
The scattering mean free path A, averaged over some 20 western events, • 
as a function of particle rigidity R. The data point at R = 1 MV refers 
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to 0. 5-1.1 MeV electrons, while those above 100 MV refer to protons 
above 4 MeV. The error bars indicate the standard deviation 
of the average, rather than the uncertainties in the measurements. 

For protons between 4 and 80 MeV, X decreases with' rigidity. 

Figure 10 Frequency distribution of the Gaussian injection width cx for 18 0.5- 
1.1 MeV electron, 30-50 MeV proton, and 9-19 MeV proton events. 

Figure ljl Five examples of proton events with prolonged injections each in a 
different energy range. Note that injection widths up to 8 hours 
have been found for proton events. 

Figure 12 The Gaussian injection width a as a function of the particle velocity 
g averaged over some 20 events. The data point at 8 = 0.9 refers to 

/ 

0.5- 1.1 MeV electrons, and the rest refer to protons. The error bars 
show the standard deviation of the average. The dashed curve gives 
the best fit in the form, a <= 8~ 0,55 . Long-inj ection events were not 
included in the averaging. 

Figure 13 The Gaussian injection width a versus 8 for two typical events. In 
each case, the relation cr “ 8 appears to be appropriate. Data 

at 8 = 0.9 and 1.0 refer to 0. 5-1.1 MeV and 3-12 MeV electrons, 
respectively. The rest refer to protons ranging from 4 to ^ 300 MeV 
in energy. 

Figure 14 The Gaussian injection width o versus 8 for two extreme cases, a 

decreases steeply with particle velocity for the March 29, 1970 event, 
while it shows no variation as a function 'of B for the April 22, 1971 


event , 
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